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' Abstract 

^ In the model of Arkani-Hamed, Dimopoulos, and Dvah where gravity is allowed to propagate in the 

D extra dimensions of very large size, virtual graviton exchange between the standard model particles can 

give rise to signatures that can be tested in collider experiments. We study these effects in dilepton and 
diphoton production at hadron colliders. Specifically, we examine the double differential cross-section 
" ^ , in the invariant mass and scattering angle, which is found to be useful in separating the gravity effects 

00 . from the standard model. In this work, sensitivity obtained using the double differential cross-section 

■ is higher than that in previous studies based on single differential distributions. Assuming no excess of 

. events over the standard model predictions, we obtain the following 95% confidence level lower limits on 

Q> . the effective Planck scale: 0.9 - 1.5 TeV in the Tevatron Run I, 1.3 - 2.5 TeV in Run Ila, 1.7 - 3.5 TeV 



in Run lib, and 6.5 — 12.8 TeV at the LHC. The range of numbers corresponds to the number of extra 
dimensions n = 7 — 2. 



:a 

^ ■ 1. Introduction 

, Recent advances in string theory have revolutionized particle phenomenology. Namely, the previously un- 

■ reachable Planck, string, and grand unification scales (Afpi, M^t, and A/gut, respectively) can now be 

brought down to a TeV range If this is the case, one expects low energy phenomenology that can be 
tested in current and future collider experiments. 

An attractive realization of the above idea was recently proposed by Arkani-Hamed, Dimopoulos, and 
Dvali In their model, the standard model (SM) particles live on a D3-brane, predicted in the string 
theory, and the SM gauge interactions are confined to this brane. On the other hand, gravity is allowed to 
propagate in the extra dimensions. In order to bring the Planck scale (10^^ GeV) to the TeV range, the 
size of these compactified dimensions is made very large compared to (Afwcak)^^- The relation among the 
Planck scale Mpi, size R of the extra dimensions, and the effective Planck scale Ms is given by: 

Af|i - Af^'+2 i?" , (1) 

where n is the number of extra (compactified) dimensions. From this relation, the size R of the compactified 
extra dimensions can be estimated. Assuming that the effective Planck scale Ms is in the TeV range, 
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it gives a very large R of the size of our solar system for n = 1, which is obviously ruled out by the 
experiment. However, for all n > 2 the expected R is less than 1 mm, and therefore does not contradict 
existing gravitational experiments. 

With the SM particles residing on the brane and the graviton freely propagating in the extra dimensions, 
the SM particles can couple to a graviton with a strength comparable to that of the electroweak interactions. 
A graviton in the extra dimensions is equivalent, from the 4D-point of view, to a tower of infinite number of 
Kaluza-Klein (KK) states with masses Mk — 2Trk/R [k = 0, 1, 2, oo). The coupling to each of these KK 
states is ~ 1/Afpi. The overall coupling is, however, obtained by summing over all the KK! states, and thus is 
^ l/Ms- Since Ms is in the TeV range, the gravitational interaction is as strong as electroweak interactions, 
and thus can give rise to many consequences that can be tested in both the accelerator and non-accelerator 
experiments. 

A large number of phenomenological studies in this area have recently appeared. Among these studies, 
the strongest lower bound on the effective Planck scale (30-100 TeV for n — 2) comes from astrophysical 
(SN1987A) and cosmological constraints Collider signals and constraints [§[ |[ |[ 0, ||, ||] come from 
diboson, dilepton, dijet, top-pair production, and real graviton emissions. 

In general, present collider experiments are sensitive to the effective Planck scale below 1 TeV. In Refs. 

the Drell-Yan and diphoton production at the Tevatron were used to constrain the scale Als- In these 
studies, however, only the invariant mass distribution of the Icpton or photon pair is used. We found that 
the distribution in the central scattering angle, in addition to the invariant mass distribution, further helps 
to constrain the scale Mg. 

In this work, we use the double differential cross-section, (Pa/dMdcos9* , to probe the effective Planck 
scale Ms in Run I and Run II at the Tevatron and at the LHC. The advantage of using double differential 
distribution is that the differences in the invariant mass and scattering angle between the SM and the 
gravity model can be contrasted simultaneously. Furthermore, for a 2 — > 2 process the invariant mass M 
and the central scattering angle cos 9* already span the entire phase space. We, therefore, do not need to 
optimize the kinematic cuts or choose optimal variables {e.g., forward-backward asymmetry, charged forward- 
backward asymmetry, etc.), because all the relevant information is already contained in the (Af x cos^*) 
distribution. We will show that sensitivity obtained in this study has improved substantially, compared to 
previous studies, in which only single differential distributions were used. By analyzing double differential 
distributions in dilepton and diphoton production simultaneously, we are able to reach sensitivity on Ms at 
the 95% confidence level (C.L.) as high as 0.9 - 1.5 TeV in the Tevatron Run I, 1.3 - 2.5 TeV in Run Ila, 
1.7 - 3.5 TeV in Run lib, and 6.5 - 12.8 TeV at the LHC, for n = 7 - 2. This is our main result. 

The organization of the paper is as follows. In the next section, we give the cross section for dilepton 
and diphoton production in the presence of strong TeV scale gravity. In Sec. |3] we describe the procedures 
in estimating the sensitivity limits. In Sec. ^ we present our results for the Tevatron and for the LHC, and 
we conclude in Sec. b] 
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2. Drell-Yan and Diphoton Production 



An effective Lagrangian for the low scale gravity interactions between the SM particles and the graviton was 
derived by Han et al. in Ref. This effective theory is valid up to a scale of about Ms- The Drell-Yan 
production, including the contributions from the SM, gravity, and the interference terms, is given by Q: 



dMudydz 
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Here ^/s is the center-of-mass energy of the pp collision, z = cos 9* is cosine of the scattering angle in 
the parton center-of-mass frame, y is the rapidity of the lepton pair, fq/g{x) is the parton distribution 



function, and we have assumed that M| ^ s, \ 



In the above equations, s — Mf^, xi^2 — "^e^^j 
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Tzf - QfSiv? 6'w, g{i = -Q/ sin^ 9^, gl = + gij)/2, and g[ = - g^)/'^- It is implied that all 
possible qq initial states are summed over. In what follows, we substitute ry = J- /Mg for convenience and 
for use as a fit parameter. 

Similarly, we calculate the diphoton production. The double differential cross section is given by 0: 



dMjjdydz 



E 



1 



ASttsM, 



fq{xi)fq{x2) 
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(4) 



where z = cos f?* is the cosine of the scattering angle in the parton center-of-mass frame and y is the rapidity 
of the photon pair. For compatibility with the Drell-Yan channel we use the range of z in Eq. (Q) from —1 
to 1, even though the final state photons are indistinguishable from each other. We account for NLO QCD 
corrections via a if-factor (see Eqs. (^) and (||)). We use K = 1.3 in the calculations. 

For the diagrams with virtual Kaluza-Klein graviton exchange it is necessary to introduce an explicit 
upper cut-off, of the order of Mg, to keep the sum over the KK states finite. A naive argument for the 
existence of such a cut-off is that the effective theory breaks down above Ms, where detailed understanding 
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of string dynamics is required. A recent observation by Bando et al. | |lO| suggests a way around this issue 
by postulating that the brane is actuaUy "flexible," with a certain tension. When the SM particles that live 
on the brane couple to the Kaluza-Klein states of a bulk gravitational field, the brane has to "stretch" out 
in order to "catch" these Kaluza-Klein states. These stretches are actually quantum fluctuations, usually 
suppressed exponentially. The higher the n of the Kaluza-Klein state, the stronger the suppression is. From 
the above argument, the contribution of high Kaluza-Klein states is suppressed, and the arbitrary cut-off in 
the sum over the KK states becomes irrelevant. Bando et al. showed that if the brane tension is equal 
to Ms, total amplitude with such a suppression is the same as the amplitude on a non- flexible brane with a 
cut-off scale set at Ms- 

Another argument, coming from the fundamental string theory, is that the coupling constant of each 
Kaluza-Klein state (Mfe = 27rfc/i?) is, in general, not independent of fc, but exponentially suppressed |pl|: 



g{k) ~ ga{k) exp ^2^/2 ■ (5) 

Here a{k) depends on the normalization of the gauge kinetic term and c is a constant that depends on 
the fundamental theory. Hence, even though all the Kaluza-Klein levels are summed over, the sum is not 
divergent and is equivalent to the unsuppressed sum with a certain cut-off scale. 

Following these recent observations, we relaxed the assumption about the cut-off scale A when summing 
the effects of all virtual graviton propagators. In Han et al., at amplitude level, the sum of the propagators 
J2k ^/{^ ~ "^fe + i"'' truncated for nik > A, where A is chosen to be Ms. In this work, we allow A to be 
different from Ms, but in order for the effective theory to be valid, it is required that A ^ 0{Ms). We define 
a scale factor c = 0{1), such that A ~ cMs. After this modification, the corresponding change in the above 
equations is: 

In the numerical analysis, we use 77 as the fit parameter in order to reduce non-linearity of the problem. 
Once the best fit value of 77 is obtained, it is straightforward to obtain the corresponding value of AIs for 
given n and c. 



3. Procedures 

3. a Experimental Acceptance 

We use typical kinematic and geometrical acceptance of the D0 detector to estimate the sensitivity in Run 
I and Run II of the Tevatron for both dilepton and diphoton production: 

\y,\ < 1.1 or 1.5 < \y^\ < 2.5 , 
Mu > 50 GeV and px^ > 25 GeV , 
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where i — fi, or 7. The integrated luminosities used in our study are 130 pb~^, 2 fb~^, and 20 fb^-'^ for 
Run I, Ila, and lib, respectively. In addition to the acceptance losses we take into account the detector 
resolution effects, as well as the longitudinal smearing of the primary interaction vertex. 
For the LHC (^/s = 14 TeV pp collision) we use the following "typical" acceptance cuts: 

\yii)\ < 2.5 , Mu > 50 GeV , pr, > 25 GeV , (i = e, ^i, 7) . (7) 

In addition to the acceptance losses, we assume the efficiency of either dilepton or diphoton reconstruction 
and identification to be 90% for the LHC or Run II of the Tevatron, and 80% for Run I. |^ In the case of 
charged leptons, detection inefficiency comes from the requirements on consistency of a track in the central 
detector with a calorimeter energy deposition (electrons) or a track in the outer muon detector (muons). For 
photons the inefficiency primarily comes from the losses due to photon conversions in the material in front 
of the central tracker. 

3.b Monte Carlo data generation 

In order to estimate the sensitivity of collider experiments to the low scale gravity model, we need to generate 
some "realistic" data sets. To set limits on the scale Ms, we assume that the SM is correct up to the energies 
of the Tevatron or the LHC. We use the SM cross section of dilepton production (the first line in Eq. (||)) to 
generate a smooth double differential distribution in Mu and z — cos 9* . We divide the Mu x z plane into a 
grid of 20 x 20 (50 x 20) bins, with AIu from GeV to 2000 (10000) GeV and z from -1 to 1 for Tevatron 
(LHG). For each bin (i, j) of this grid, the expected number of events, SfJ^, is obtained by multiplying the 
cross section in this bin by the known integrated luminosity and efficiency. We further proceed with a Monte 
Carlo (MG) gedankenexperiment. For each bin (j,j) we generate a random number of events, riij, using 
Poisson statistics with the mean SfJ^. Similar gedankenexperiment is done for the diphoton production. We 
use the dilepton or diphoton MC data sets generated in this way to perform the best fit to the low scale 
gravity model (see section Either of the two channels, or their combination can be used in the fit. |^ 

3.C Fitting procedure 

We extract the lower limit on the gravity scale Ms by fitting the "data" obtained in a MC experiment 
with a sum of the SM background and Kaluza-Klein graviton contribution. We employ both the maximum 
likelihood method and pure Bayesian approach with a fiat prior probability for ?7 > and for 77 < 0. Since 

*Dimuon acceptance and efficiency were significantly lower in Run I, but we deliberately have not done a more realistic 

simulation in this case. First, the contribution of the dimuon channel to the overall sensitivity is very small. Second, a 

designated data analysis, which is currently being finalized by the D0 Collaboration jl2) will soon override our estimates by 

utilizing real collider data and measured efficiencies. 

tNote that combination of the dilepton and diphoton channels implies combination of the corresponding likelihoods, not the 

spectra! 
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we focus on the number of large extra dimensions > 3, the fitting procedure is straightforward, as the 
factor T can be taken out of the integration over the phase space. 

We generate three templates that describe the cross section in the case of large extra dimensions. The 
first one describes the SM cross section on the rectangular grid described in section [3.b| . The other two 
describe terms proportional to 77 (interference term) and to 77^ (Kaluza-Klein term), respectively. We then 
parameterize production cross section in each bin of the Af x z grid as a bilinear form in r\: 

O- = CTSM + Cr4'7 + 0-8??^, (8) 

where ctsm, "'4, and cts are the three templates described above. 

In Figs. 1^ and|^, we show the 3-D plots for the pure SM, the interference, and pure gravity contributions 
for dilepton and diphoton production, respectively, at the 2 TeV Tevatron. It is clear that the pure SM 
decreases rapidly with the invariant mass. This is in contrast with the pure gravity contribution that rises 
quite sharply with the invariant mass and then turns over due to the effect of parton distribution functions. 
The interference term also shows similar characteristics. The angular distribution also exhibits substantial 
difference among the pure SM, pure gravity, and the interference. Note the asymmetry of the interference 
term for dilepton production (Fig. |l|b) that arises from the charge asymmetry of the Tevatron beams and 
final state particles. Analogous distribution for diphotons or in the LHC case is symmetric. 

The probability to observe certain set of data M = {f^y}, where («,j) are the bins in M and cos0*, 
respectively, as a function of 77 is given by the Poisson statistics: 



p{m^Y^ ^''^'\ (9) 



a'' 

e 



where Sij = Ctaij, and C is the integrated luminosity, e is the identification efficiency, and Uij is the cross 
section given by Eq. (||), integrated over the bin (j,j). 

We now can use Bayes theorem to obtain the probability of 77, given the observed set M: 



= jfdx exp (^^^ ) P{J^\V), (10) 



where A is the normalization constant, obtained from the unitarity requirement: 

dv P{v\^f) - 1 , (11) 





xq is the central value of the eC, and ax is the assumed Gaussian error on the quantity x. In order to 
minimize the uncertainty we perform in situ calibration by normalizing x to reproduce the observed 
number of events with M < 100 GeV (200 GeV) at the Tevatron (LHG) {i.e., we use the first mass bin of the 
MG grid to perform the normalization). Such a procedure is justified by the fact that possible contribution 
from Kaluza-Klein gravitons virtually does not affect the low mass region (see Figs. |l| and ^). We, therefore, 

tPor the case of n = 2, depends on s. In the next section, wo will also give the results for n = 2 by estimating the average 
s for the gravity term in dilepton and diphoton production. 
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assume cr^ to be 10% or (100%/y "-ij); whichever is smaller. (When setting limits on 77 we then only use 
the mass bins above the normalization region, i.e. i > 1.) 

The 95% C.L. limit on signal, 7795, is obtained from the following integral equation: 

/ driP{r]\Af) = 0.95. (12) 
Jo 

A less sophisticated likelihood approach does ignore systematic error on the efficiency and integrated 
luminosity and simply treats P{Af\ri) as the likelihood function. The 95% C.L. limit in this case is obtained 
by requiring the integral of the likelihood function from the physics boundary (77 — 0) to 7795 to be equal to 
0.95. As was mentioned before, both approaches yield very close limits on 77. While the Bayesian technique is 
a natural way to account for the systematic errors on the efficiency (and background) estimates (and this is 
the approach actually used by the D0 experiment to derive limits), we implemented the classical likelihood 
approach as well, primarily to demonstrate the robustness of the limit setting technique. 

We further combine the results obtained from the dilepton and diphoton channels by adding the prob- 
abilities (likelihoods) and solving the integral equation (^2|) (or its equivalent for the maximum likelihood 
method). 

As an additional cross check we have tested the fitting techniques with a set of the MC experiments 
assuming a non-zero Kaluza-Klein graviton contribution. Both the Bayesian and maximum likelihood fits 
were capable of extracting the input value of the gravity scale without a systematic bias, as expected. 

To convert 7/95 from a single MC experiment into a measure of sensitivity of future experiments, we 
repeat the above procedures (both the gedankenexperiment and fit) many times. The limits obtained in 
these repeating experiments are histogrammed. Sensitivity to the parameter 77 is defined as the median of 
this histogram, i.e. the point on the sensitivity curve which 50% of future experiments will exceed. All the 
limits given in the next section are based on this sensitivity measure. (An alternative approach that defines 
sensitivity as the most probable outcome of the gedankenexperiment agrees with the one we used within 5% 
accuracy.) 

4. Results 

In our study, we include both the electron and muon channels in the Drell-Yan production. In Table |l| we 
show the sensitivity to 77 in Run I, Run II of the Tevatron, and at the LHC using dilepton and diphoton 
production, as well as their combination. Corresponding Mg reach is also shown for n ~ 2 — 7 and c = 1. 
For other values of c they are shown in Fig. ^ or can be calculated by simple rescaling, using Eq. (^ . For 
the case n = 2 the conversion of 77 limits into A/5 limits is not straight- forward, as it depends on the s of the 
subprocess, see Eq. We use the pure gravity contribution in the dilepton and diphoton production to 
estimate the corresponding average s. With the average s we can then roughly estimate the Ms limits for 
n ^ 2. For diphoton production the average s for Run I, Run II, and LHC are (0.61 TeV)^, (0.66 TeV)^, 
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Table 1: Sensitivity to the low scale gravity model parameter 77 — ^c"~^/M| in Run I, Run II of the Tevatron and 
at the LHC, using the dilepton, diphoton production, and their combination. The corresponding 95% C.L. limits on 
Ms are given in TeV for n — 2 — 7 and c = 1. Results for other c values can be read from Fig. ^ or obtained by 
rescaling, using Eq. (^. 





7/95 (TcV-4) 


n = 2 


n = 3 


n = 4 


n = 5 


71 = 6 


n = 7 


Run I (130 pb-i) 
















Dilepton 


0.66 


1.21 


1.32 


1.11 


1.00 


0.93 


0.88 


Diphoton 


0.44 


1.39 


1.46 


1.23 


1.11 


1.03 


0.98 


Combined 


0.37 


1.48 


1.53 


1.29 


1.16 


1.08 


1.02 


Run Ila (2 fb^i) 
















Dilepton 


0.163 


1.92 


1.87 


1.57 


1.42 


1.32 


1.25 


Diphoton 


0.077 


2.40 


2.26 


1.90 


1.71 


1.60 


1.51 


Combined 


0.072 


2.46 


2.30 


1.93 


1.74 


1.62 


1.54 


Run lib (20 fb-i) 
















Dilepton 


0.054 


2.70 


2.47 


2.08 


1.88 


1.75 


1.65 


Diphoton 


0.025 


3.40 


3.00 


2.53 


2.28 


2.12 


2.01 


Combined 


0.021 


3.54 


3.11 


2.61 


2.36 


2.20 


2.08 


LHC (14 TeV, 100 fb^i) 
















Dilepton 


2.20 X 10-* 


10.2 


9.76 


8.21 


7.42 


6.90 


6.53 


Diphoton 


1.24 X lO-'' 


12.1 


11.3 


9.47 


8.56 


7.97 


7.53 


Combined 


1.05 X 10"" 


12.8 


11.7 


9.87 


8.92 


8.30 


7.85 



and (3.2 TeV)^, respectively, while for dilepton production the average s are (0.60 TeV)^, (0.64 TeV)^, and 
(3.1 TeV)2, respectively. 

The Drell-Yan channel is not as sensitive as the diphoton channel and, therefore, the combined limit is 
close to the limit from the diphoton channel only. In Run I, using the combination of two channels, the 
sensitivity to Mg is about 1.0 to 1.5 TeV for ti = 7 — 2 and c — 1. It increases to 1.5 to 2.5 TeV in Run 
Ila, and 2.1 to 3.5 TeV in Run lib. At the LHC, the sensitivity soars up to 7.9 - 12.8 TeV. Both higher 
center-of-mass energy and increase in the integrated luminosity help to improve the limits. 

We also study the improvement in the sensitivity from the double differential (Pa/dMd cos 6* fit compared 
to that from the single differential da/dM fit. We have repeated the entire procedure with a 20 x 1 grid in the 
{M X COS0*) plane, which is equivalent to fitting the single differential distribution da/dM. Corresponding 
limits in Run Ila deteriorate to: 

7795 = 0.176 TeV"'' for dilepton , (13) 
T]95 = 0.089 TeV* for diphoton , (14) 
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^795 = 0.084 TeV for combined dilepton and diphoton . (15) 

By using the double differential cross section we achieve an improvement of about 10% (15%) in the limit 
on rj for dileptons (diphotons) . While such an improvement in sensitivity translates only into a few per cent 
increase in the limit on AIs, it is actually equivalent to a 30% decrease in the integrated luminosity, required 
to set a certain limit on Ms- 

5. Conclusions 

The sensitivity to the effective Planck scale Ms obtained in this analysis supercedes those from the previ- 
ous studies, in which only one-dimensional distributions were used {e.g., Drell-Yan production diboson 
production diphoton production |j, dijet production and top pair production j^). The recent work 
by Eboli et al. that studied diphoton production in the Tevatron Run Ila and at the LHC quotes 95% 
C.L. upper hmits on Ms of 1.73 TeV (n = 4) in Run Ila and 7.7 TeV {n = 4) at the LHC. Our hmits 
exceeds the latter, partly because we have taken into account the invariant mass and angular distributions 
simultaneously, and partly because we do not impose the unitarity constraint -x/I < 0.9Ms and use a slightly 
higher efficiency. 

As we have mentioned in the Introduction, the invariant mass M and the central scattering angle cos 6* 
already span the entire phase space of a 2 — + 2 process. Thus, our fit method gives an ultimate way of 
probing the low scale gravity in the virtual graviton exchange processes, because all relevant information 
is contained in the {M x cob 6*) plane. We have shown that the improvement in the limits of rj from the 
double differential d? a / dM d cos 9* fit over those from the single differential da/dM fit is about 15%, which 
corresponds to a 30% decrease in the integrated luminosity needed to obtain a certain sensitivity in Ms- 

To summarize, we have analyzed the double differential distribution in the invariant mass and scattering 
angle for dilepton and diphoton production at hadron colliders. We have obtained better sensitivity than 
previous studies have achieved. Limits that we obtained using the Bayesian approach and maximum like- 
lihood method are numerically identical. The expected limits on Ms are: 0.9 — 1.5 TeV (Run I), 1.3 — 2.5 
TeV (Run Ila), 1.7 - 3.5 TeV (Run lib), and 6.5 - 12.8 TeV (LHC) for n = 7 - 2. 
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Figure 1: The 3-D plots for double differential distribution d^a /dMud cos 9* for Drell-Yan production at the 2 TeV 
Tevatron. (a) SM only, (b) the interference term between the SM and the gravity contributions, proportional to 
and (c) the pure gravity contribution, proportional to . Here r) = J^c"~^ /Mg. Note that in (b) linear scale is used 
in the 0-axis in order to show the negative z values. 
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Figure 2: The 3-D plots for double differential distribution a /dM-y-yd cos 0* for diphoton production at the 2 TeV 
Tevatron. (a) SM only, (b) the interference terms between the SM and the gravity contributions, proportional to 77, 
and (c) the pure gravity contribution, proportional to rf. Here rj = J^c"~^ /Mg. 



13 



10' 



10" 



> 

CD 



10 ■ 



10 ■ 



10' 



> 

CD 



10"- 



10 ■ 



10' 



■p 10" 
> 



10 ■ 



10"' 




Ms (TeV) 



Figure 3: rj versus Ms for a given (n, c). For c — 0.5, 1 the lines from top to bottom are for n = 3, 5, 7, whereas for 

c = 2 the Unes from top to bottom are for n — 7,5, 3. The 95% C.L. hmit on 77 for Tevatron Run I, Run II, and for 

LHC are shown. Limits are based on: (a) Drell-Yan production, (b) diphoton production, and (c) combined. 
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